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Thermodynamic behavior and stability of Polytropic gas 
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We focus on the thermodynamic behavior of Polytropic gas as a candidate for dark energy. We use 
the general arguments of thermodynamics to investigate its properties and behavior. We find that a 
Polytropic gas may exhibit the dark energy like behavior in the large volume and low temperature 
limits. It also may be used to simulate a fluid with zero pressure at the small volume and high 
temperature limits. Briefly, our study shows that this gas may be used to describe the universe 
expansion history from the matter dominated era to the current accelerating era. By applying some 
initial condition to the system, we can establish a relation between the Polytropic gas parameters 
and initial conditions. Relationships with related works has also been addressed. 


I. INTRODUCTION 

In the Einstein General Relativity framework, 
Friedmann-Robertson-Walker metric (FRW) and its con¬ 
formal form help us get a suitable model for describ¬ 
ing the universe expansion and its inhomogeneities in 
scales smaller than 100—Mpc, respectively [JQ- In stan¬ 
dard cosmology, the universe is born from a singular¬ 
ity called big bang and is inflated in primary moments. 
Thereinafter, the rate of this expansion is decreased, 
since Radiation and Matter play the role of dominated 
fluid, which is the determinant of the universe expan¬ 
sion rate, in subsequent eras [l|. Finally, the universe 
deals with an accelerated expansion with positive rate, 
about 13700—Myr after big bang [3, [3-[l|- Although the 
latest and current phase of universe expansion satisfies 
the thermodynamics stability conditions Ml but, the 
nature of dominated fluid, which supports this phase of 
expansion and called dark energy (DE), is a mysterious 
puzzle [iMl. 

In addition to the mentioned puzzle, standard cosmol¬ 
ogy suffers from a series of other weaknesses such as 
the coincidence and fine tuning problems [l|. There are 
various attempts to give a candidate for describing DE 
and solving the standard cosmology problems |12 - [l4l |. 
Among these attempts, some models include a varying 
DE candidate with equation of state P = p + f{p), where 
P and p are the DE pressure and energy density, respec¬ 
tively [Illll. In astrophysics, a Polytropic gas with 
equation of state 


has numerous applications [I3- In addition, n called the 
Polytropic index and AT is a constant iTp. An adiabatic 
gas [l^ and a degenerate electron gas [l7j| are two ex¬ 
amples of such gases. Since the equation of state of this 
gas is in the P = f{p) form, it attracts the attention of 
cosmologist to itself as a probable way to describe DE 
and the coincidence as well as the fine tuning problems 


Karami et al. investigated the mutual interaction be¬ 
tween Polytropic gas (as the DE candidate) and CDM, 
and find out that for some even and positive Polytropic 
index and some positive values of K, the Poly tropic gas 
behaves as the Phantom dark energy [l^. Moreover, the 
generalized second law of thermodynamics is always sat¬ 
isfied by a universe filled with a Polytropic gas model of 
DE together with a CDM which interact with each other 
. Authors in [2l| used the first law of thermodynam¬ 
ics, and qualitatively show that a Polytropic gas may 
exhibit a DE behavior. In addition, they show that this 
model leads to a suitable fitting with observational data 
about the current expanding era [2l|. In fact, Asadzadeh 
et al. have argued that a Polytropic gas model of DE 
with Polytropic index n < — 1 and A' < 0 is in suitable 
agreement with the observational data 22|. It is also 
shown that the Polytropic gas model of DE may simu¬ 
late the new agegraphic as well as the holographic dark 
energy models and thus the current accelerating phase of 
universe [1^ . Interaction between the Polytropic gas 

model of DE and dark matter in Kaluza-Klein cosmology 
has also been studied 1^. 


P = Kp^^+i\ ( 1 ) 

where 


is its energy density, while U and V are the total en¬ 
ergy of gas and the volume of its container, respectively, 
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It is useful to mention here that in all of Refs. [H- 
1^ , authors investigate the Polytropic gas properties and 
behaviors by considering the cosmological setup. As a 
common result, all of them find that, in the FRW back¬ 
ground, the Polytropic gas may be a suitable candidate 
for describing DE. Therefore, one may ask that what is 
the origin of the DE like behavior of Polytropic gas? Does 
this ability of the Polytropic gas to describe DE is due 
to the gravitational theory used to describe the universe 
expansion? Indeed, does thermodynamical properties of 
Polytropic gas is similar with those of DE without consid¬ 
ering a cosmological background? Such confusions have 
been arisen in studying the Chaplygin gas as well as its 
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generalized and modified forms [11,13^. In fact, authors 
in by focusing on the thermodynamic arguments, 

have been shown that such gases have enough ability for 
exhibiting the DE like behavior. Although, the equation 
of state of Polytropic gas is similar to that of the gen¬ 
eralized Chaplygin gas but the physics behind them are 
completely different. Whiles, the Chaplygin gas model 
comes from the Russian physicist attempts to study the 
lifting force which applies on the plane wings in aerody¬ 
namics [13 , the Polytropic gas models are introduced to 
study the hydrostatic equilibrium equation and thus the 
star evolution [i3- Therefore, it is not straightforward 
to generalize the results of Refs. to Polytropic gas 

models. 

Here, we want to study the thermodynamics of a Poly¬ 
tropic gas confined to volume V. In fact, we are eager 
to know that under which thermodynamic conditions a 
Polytropic gas exhibits the DE like behavior. Does this 
model satisfies the thermodynamic stability conditions, 
and does it pleases the thermodynamic expectations? 
Our study shows that the Polytropic gas may naturally 
exhibit the DE like behavior. Since we only use thermo¬ 
dynamics to study the behavior of Polytropic gas, our 
results are independent of considering the gravitational 
theory used to describe the universe expansion. 

The paper is organized as follows. In the next section, 
we consider an adiabatic Polytropic gas and evaluate its 
thermodynamical properties such as its energy as a func¬ 
tion of volume. The thermodynamic stability conditions 
are also addressed. Section III includes more debates 
about the thermodynamic stability conditions and some 
examples. In section IV, by applying initial conditions 
to the system we get a relation between the Polytropic 
index, K and initial conditions. Section V includes sum¬ 
mary and concluding remarks. 


II. ADIABATIC POLYTROPIC GAS EOS 

The pressure of a fluid with energy U confined into a 
cylinder with volume V is evaluated as [l^ 


and simple calculations, we get 

D = (-l)-"iW"E(l + (A)-^)-, (6) 


which leads to 


p = ^ = (-irA-"(i + (^)^)-", (7) 

for the energy density of Poly tropic gas. This equation 
leads to 


p^{-irK-\ ( 8 ) 

and 

(9) 

for n > 0 and n < 0, respectively, in the small volume 
limits. Since it seems that density should be positive for 
even n [Hill,!!!. these equations state that 6 should 
meet the 5 > 0 condition. Since ^ > 0, one can use 
Eq. d?]) to conclude that, independent of n, density is 
positive for K < 0. Moreover, inserting this equation 
into o to reach 

P = + (10) 

0 

which yields 


P = 




and 


P 1 


( 11 ) 


( 12 ) 


as the corresponding pressure and the state parameter of 
Polytropic gas, respectively. Therefore, the state param¬ 
eter of Polytropic gas meets the — 1 < w < 0 condition. 
It is interesting to note that for b = 0, one obtains 


Vw 


s 


-P. 


( 3 ) 


Combining Eqs. © and and inserting the result 
into to get 



By taking integral from this equation, one finds the en¬ 
ergy of Polytropic gas as 

t7=(-l)-"(AE-i+6)"”, (5) 

where b = b{S). It is also apparent that Eqs. dH) and dS]) 
are also available, ii K = K{S). By defining 6~^ = -^ 


P = -p= (-1)"+^A-", 

(13) 

P 

(14) 

UJ = — = —1, 


P 


which are independent of the system volume. Therefore, 
the b = 0 case is capable to explain the primary infla¬ 
tionary era and the current expanding phase as well as 
the anti de-Sitter spacetime. For example, consider a sit¬ 
uation in which K > 0. While ^ = 0, on one hand, 
for even n this equation is compatible with the equation 
of state of the cosmological constant (p > 0) which may 
support the de-Sitter spacetime. On the other hand, an 
odd n points to a fluid with p < 0 which may support 
the anti de-Sitter spacetime. 
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Whenever n > 0 and S does not diverge, one gets 


c^-p, (15) 


for small volumes, which is the same as the result of the 
6 = 0 case uni). In order to derive this equation, we used 
Eqs. © and m- This result is in accordance with the 
primary inflationary era. It is interesting to note that the 
generalized Chaplygin gas has similar behavior at large 
volumes [s^. For large volumes, Eqs. o and (HU lead 
to 


P 


V 


and 


P 


(_l)n+l/C-n^ 




(16) 


(17) 


small volumes) to w ~ — 1 (for large volumes), it may be 
useful to describe the universe expansion from the matter 
dominated era to the current expanding phase. 

The thermodynamic stability condition or the convex¬ 
ity of the energy surface requires that (iU| 


,dP, 

(^)5<0, 


( 21 ) 


and 


Cp>Cv> 0, (22) 


where Cy and Cp are the heat capacities at constant 
volume and pressure, respectively. Additionally, since 
for a fluid with N particle 


Cp-Cv 


TVa^ 
Nkt ’ 


(23) 


respectively. We should note that since V ^ I and n > 
0, w ~ 0. In fact, while length is proportional to the 
cosmological scale factor a, V oc and we get 


P 


{-lYK- 


(18) 


where kt = — Ai, and a are the isothermal com- 
V{gv>T 

pressibility and the coefficient of thermal expansion, re¬ 
spectively [3, a Kt > 0 and Cy > 0 are simultaneously 
satisfied condition (l22l) is also met. The latter means a 
thermodynamic system, which satisfies Eq. (EU, 


which indicates a pressureless matter. This behavior of 
Polytropic gas at large volume is the same as that of the 
generalized Chaplygin gas in small volumes [s^. There¬ 
fore, a Polytropic gas with n > 0 and 6 yf 0 cannot sim¬ 
ulate a fluid with P = —p 0 at large volumes. Briefly, 
a Polytropic gas with n > 0, while 5 does not diverge, 
behaves as a fluid with a; ~ — 1 and a pressureless matter 
at the small and large volumes limits, respectively. 

Meanwhile, while n < 0 and <5 does not diverge, by 
using Eqs. 0 and (ED one gets 


Cy > 0, 


and 


dP 


(24) 


(25) 


also meets the thermodynamic stability condition. For 
the 6 = 0 case, from Eq. ED, simple calculations lead to 
(|T)g = 0 and therefore, the stability condition (1^ is 
marginally satisfied. Now, using Eq. ED to reach 


and 


P^O, p 


(-i)”7s:-" 


P Pi ~ -p, 


(19) 


( 20 ) 


as the pressure and density of Polytropic gas in small 
and large volumes, respectively. In order to get ED , we 
assumed again P oc a^. Therefore, a Poly tropic gas with 
n < 0 behaves as the pressureless matter and a fluid with 
a; ~ — 1 at small and large volumes limits, respectively. 
This behavior is fully consistent with that of the gener¬ 
alized Chaplygin gas [35l| . Let us summarize the above 
results. The 6 = 0 case indicates a fluid with oj = —1 
and therefore, may be used to simulate the primary and 
current inflationary eras. Moreover, since the state pa¬ 
rameter of a Polytropic gas with n > 0 is increased from 
OJ ~ — 1 (for small volumes) to oj ~ 0 (for large volumes), 
it may be useful to describe an expanding universe which 
expands from an inflationary like primary era to a matter 
dominated like era. Finally, since the state parameter of 
a Polytropic gas with n < 0 is decreased from w ~ 0 (for 


l + {5V)-i ■ 


(26) 


Bearing the 6 > 0 condition together with Eqs. o 
and ED in mind, a Polytropic gas with p > 0 satisfies 
condition (ED, if — 1 < n < 0. In addition, a Polytropic 
gas with negative density (p < 0) satisfies condition (1^ 
if its Polytropic index meets either the n<—lorn>0 
conditions. It is also apparent that, just the same as the 
6 = 0 case, the n = — 1 case marginally satisfies condi¬ 
tion ED- In the next section, we investigate the quality 
of validity of Eqs. (ED 


III. THERMAL POLYTROPIC GAS EOS 

One should determine the thermal equation of state 
P = P[T,V) to investigate the quality of validity of 
Eq. (1^ . In addition, since 


Cv = T{^)y, 


(27) 
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we need also to determine the thermal equation of state 
S = S{T, V) in order to study the behavior of Cy [l^ - 
Here, since in cosmolo gica l setups authors generally set 
i^T to a constant value [13, [la, j consider the 

situation in which ^ = 0 which means that K is con¬ 
stant. In order to evaluate the temperature of Polytropic 
gas, inserting Eq. (0) into 

dU 

T=(^)y, (28) 

and get 

T = {-lY+^nV^+^{K (29) 

do 

as the temperature of Polytropic gas. Bearing Eqs. © 
and (HU together with the definition of 6 in mind, this 
equation can be written as 


T=-n 


pV^+i db 

(l + (f)-)^' 


(30) 


In order to continue our analysis, we need an exact form 
for b. To achieve this aim, bearing Eq. m in mind, a 
dimensional analysis shows 


[b]-- = [U]. (31) 

In addition, from Eq. (sni) we have db oc dS. Finally, since 
\U\ = [TS] in thermodynamics [l^ . a primary simple 
election which satisfies Eq. m is 


b=ins)-i, (32) 


as the last thermal equation of state (P = P{T, V)) 
needed to investigate the system. From this equation it is 
apparent that (f^)T = 0 which indicates condition (l25ll 
is marginally satisfied by the Polytropic gas. Finally, 
since P = PiT), (|^)t and (|^)t are also zero which 
means that there is no critical point in this situation [li¬ 
lt is useful to note that this conclusion is the direct result 
of using Eq. (l32l) to obtain b. By inserting Eq. dSll) into 
Eqs. © and ©, we get 


and 


P = 


U = 


i-K)- 


[l + ((^)^(-l)^-l)-l]r. 

_ {-K)-'^V _ 

[l + ((^)^(-l)XTT-l)-l]. 


(38) 


(39) 


as the energy density and total energy of Polytropic gas, 
respectively. As we know, the entropy of a thermodynam¬ 
ical system should be positive [T^. Moreover, Eqs. dSU) 
and (1381) can be used to get 


oj = -- 


l + ((^)—(-1)—-1)- 


(40) 


as another relation for the state parameter of Polytropic 
gas. This equation indicates that for a Poly tropic gas 
with Poly tropic index satisfying the (—1)'^ = 1 condi¬ 
tion and temperature T, which meets the 0 < T < T, 
range, the state parameter encounters the — 1 < w < 0 
range. By inserting this equation into Eq. dSH), one reach 


where P* is a universal constant with temperature di¬ 
mension, which should be evaluated from other parts 
of physics such as statistical mechanics or experimental 
data. Such analysis yields similar results for generalized 
Chaplygin gas and can be found in Ref. [S^- Taking 
derivative with respect to S from this equation to get 


db 

dS 



n 


s- 


1 


Inserting this result into (pOl) to obtain 


T = 


pV^+k 

(l + (f)") 


p; 



Now, inserting Eqs. © and dSU into (IMl) to get 


(33) 


(34) 


P= (-l)"pi+i(P, "S'-i-i)[i^-hP* 

(35) 

which leads to 

(36) 


for the entropy of Poly tropic gas. By combining 
Eps. (|1 ip . (|34l) and. dMD one gets 


P = if-"(-l)"+i(l + - l]-i)-("+i), 

(37) 


S = 


V 


(— —y 


(41) 


for entropy. Hence, since from Eq. (HU w meets the — 1 < 
w < 0 condition, the S' > 0 condition is met if X^P* > 0. 
Now, using Eq. to get 


Cy = n 


;^(-i)^(E)yT[(E)yT(_i)^_i]-i, 


(42) 

for the heat capacity at constant volume, which can be 
rewritten as 


Cy = 


y(i + (^)b^ 

^ 1 

Kpr^ + l 


(43) 


where we have used Eqs. dal and dMl) to obtain this 
equation. Now, inserting Eq. m into this equation to 
obtain 


Cy = n 


ni + (T)-)^ 




w ■ 


(44) 


Therefore, the Cy > 0 condition is satisfied if 

- - — —|— > 0. By combining this result with the re- 

suit obtained from the S > 0 condition, we get if the 
—> 0 condition is satisfied, then the S > 0 and 

K p ^+1 

Cy > 0 conditions are simultaneously satisfied. 
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Some examples 


Consider a Polytropic gas with n = — , where k 

is an integer and positive number, which leads to p > 0 . 
From Eq. (l26l) . it is obvious that this gas satisfies con¬ 
dition (| 2 T|l . Eq. (1^ indicates T* > 0 to preserve the 
S' > 0 condition which yields & > 0 (l32)l . In addition, 
from Eq. , we see that the Cy > 0 condition leads to 
K < 0. Bearing Eq. (|40l) in mind, since Eq. (IT^ indicates 
— 1 < uj < 0 , we conclude that temperature should meet 
the 0 < T < T* condition. Therefore, this Polytropic gas 
behaves as a pressureless matter for T ^ T, and small 
volumes. Moreover, its state parameter and temperature 
are decreased to —1 and 0 , respectively, by increasing the 
volume. During this process, the Polytropic gas satisfies 
the thermodynamic stability condition. Loosely speak¬ 
ing, this behavior is similar to the universe expansion 
history from the matter dominated era to the current 
accelerating phase [l|. It is useful to note that the gen¬ 
eralized Chaplygin ga,s has also similar behavior as that 
of mentioned here 13511. 


As another example, focus on a Polytropic gas with 
positive density, which yields P < 0 , and n = — 
where k is again a positive integer number. Since we 
assumed p > 0, Eq. 0 implies AT < 0. Eq. ([26|) shows 
that condition ((^T]) is fulfilled. Bearing the S definition in 
mind, in order to get positive values for S, b should meet 
the b < 0 condition. Considering Eq. (l32|) . the latter 
leads to T* < 0. From Eqs. (HU) and dm) it is apparent 
that the S > 0 and Cy > 0 conditions are simultaneously 
met in this situation, respectively. Bearing Eq. dm in 
mind, in order to preserve the — 1 < w < 0 condition, 
which comes from Eq. dm, one gets 0 < T < —T* as 
a permissible range for temperature. The latter is fully 
consistent with physical expectance about positivity of 
T. Finally, we may conclude that the investigated Poly¬ 
tropic gas may provide a suitable model for describing the 
sources of universe expansion in the matter and current 
accelerated eras. 


Finally, consider a Polytropic gas with n < — 1 and 
AT < 0. According to 0 and dm, it is obvious that 
p > 0 and P < 0, respectively. From Eq. dm, it is ap¬ 
parent that condition dm is not satisfied. Also, it is easy 
to show that the S > 0 and Cy > 0 conditions may be 
met. For example, if T* > 0 and AT” > 0 then the S' > 0 
and Cy > 0 conditions are satisfied. Bearing Eqs. dm 
and (im ill mind, this gas behaves as a pressureless mat¬ 
ter and a fluid with a; ~ — 1 for small volumes and at the 
T ^ Ttt limit and for large volumes and at the T ^ 0 
limit, respectively. It is shown that such a Poly tropic gas 
may lead to a satisfactory fitting with the observational 
data (m| . 


IV. THE POLYTROPIC GAS PARAMETERS 
FROM THERMODYNAMIC ARGUMENTS 


Consider a Polytropic gas with initial condition V = 
Vb, P = Po, p = Po and T = Pq- From Eq. dS]) we get 

^=-(Po”" (45) 

Now, inserting this result into Eqs. 0 and o to obtain 

p = Po(-Apo"+(l + Apo")(^)i)-", (46) 

ro 

and 

P = Api+^ = ApJ+"(-Apo^ + (l+Apo^)(-^)^)-(-+i), 

ro 

(47) 

for the density and pressure of Poly tropic gas, respec¬ 
tively. Here, we define new parameters e, u, p, 7 , t and 
p as 


e 

7 


—, 

Po 

-KpS , 


V _ P 
Vo’ a :-"’ 




(48) 


It is useful to mention that 7 is nothing but K defined 
in ( 2 ^. Therefore, Eqs. (j46p cinci (I3TI) Ccin. be written ns 

e = (7 +(1- 7 )n")"”, (49) 


and 

P = [(-(1 + [(|)^(-1)^ - l])-i)]-("+i), (50) 

respectively. Moreover, from Eq. dm we obtain 


P = 


P 


= K 


n+1 


n + 1 {^-lY+l^n + l 

p " = 


(7 -|- (1 — 7 )Un )"'+! 


(51) 


Inserting initial conditions into Eqs. dm and dm, and 
equating them to each other to get 


- 


(- 1 )^ 


(I - 7 )"+!' 


(52) 


In fact, this equation helps us to relate AT, n, T*, To and 
Po to each other. Now, consider a Polytropic gas with 
p > 0 and n = — , where k is an integer and positive 

number. Now, combining Eqs. (HSl) and (15^ to obtain 


K = [&^ - l]pl'"^. 
1 * 


(53) 


Therefore, for 0 < Tq < T*, since A: > 0, this equation in¬ 
dicates K < 0 which is in line with the results obtained in 
the previous section. Indeed, for every Polytropic index 
n we get 


K = 


(f)^(-I)^-I 

Po 


(54) 
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which indicates that in order to get real values for K, 
whiles T* > 0 , n should meet the (— 1 )^^ = ±1 con¬ 
dition. For example, if 0 < Tq < T* and n meets the 
(— 1 )— = 1 condition, we get K < 0 and thus 7 > 0 
meaning that K > 0 which is fully consistent with the 
results obtained in Ref. [l^. Loosely speaking, by apply¬ 
ing some initial conditions on the Polytropic gas we get 
a relation between the initial conditions and the Poly¬ 
tropic gas parameters including the Polytropic index (n) 
and constant K. 


V. SUMMARY AND CONCLUDING REMARKS 

Here, we have considered an adiabatic Polytropic gas 
and found out its energy density and pressure by using 
pressure definition in thermodynamics and the equation 
of state of Polytropic gas, respectively. Thereinafter, we 
got the corresponding relations for the state parameter 
as well as the total energy of Polytropic gas. In addi¬ 
tion, by investigating the asymptotic behavior of energy 
density of Poly tropic gas for small and large volumes, 
and using the positivit y of its density for even Polytropic 
indexes condition [ 13 , [l^ [ZJ , we got the (5 > 0 condi¬ 
tion for parameter S appeared in relations. Our study 
shows that the state parameter of Polytropic gas meets 
the — 1 < w < 0 range for the finite values of 6. In 
continue, we reviewed the thermodynamic stability con¬ 
dition which comes from the energy surfaces convexity 
condition [l^. In order to investigate all of the stability 
conditions, we used dimensional analysis to get a pri¬ 
mary relation for b and thus temperature. Due to our 
special choice for temperature, no critical point observed 
for Polytropic gas. We should note that this result may 
be changed by choosing different functions for b. We 
also pointed to the entropy of Polytropic gas, and get to 
the > 0 condition to have positive entropy. Our 


study shows that whenever the - - —j— > 0 condi- 

tion is satisfied, the Cy > 0 condition is also met by 
Polytropic gas. We have pointed to the behavior of some 
Poly tropic gases, and found out that a Poly tropic gas 
with n = — 2 k+i > ’"^here k is a. positive integer number, 
and p> 0 satisfies all of the stability conditions if T* >0 
and K < 0. Under such criterion, the state parameter 
and temperature of Polytropic gas meet the — 1 < w < 0 
and 0 < T <T^, ranges, respectively. Polytropic gas with 
n = — 2 k +3 ^ '"^tiere A: is a positive integer number, and 
AT < 0 is also addressed. We found out that the S' > 0 
and Cv > 0 conditions indicate T, < 0. Finally, we saw 
that its temperature should meet the 0 < T < —T, range 
in order to preserve the — 1 < w < 0 condition obtained 
from Eq. Our study also shows that a Polytropic 

gas with n < — 1 and K < 0 may satisfy the S > 0 and 
Cv > 0 conditions whiles T* > 0 and S is considered as 
a finite and well-defined quantity. Such a Polytropic gas 
has positive density and behaves as a pressureless mat¬ 
ter for small volumes and at the T —^ T, limit and a 
fluid with w —>■ —1 for large volumes and at the T —)> 0 
limit, whiles A" >0. It is worthwhile to mention here 
that a Polytropic gas with n < — 1 and K < 0 has an 
appropriate fitting with observational data p^ . Finally, 
by imposing some initial conditions to under investiga¬ 
tion Polytropic gas, we could find a relation between the 
Polytropic gas parameters, including n and K, and the 
initial conditions applied on system. 
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